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Abstract 

Transonic aeroelastic stability analyses and flutter 
calculations are presented for a generic Uansport-type wing 
based on the use of the CAP-TSD (Computational 
A.eroelasticity Rrogram - Iransonic Rmall Qjsturbance) 
finite-difference code. The CAP-TSD code has been recently 
developed for transonic unsteady aerodynamic and aeroelastic 
analysis of complete aircraft configurations. In this study, a 
binary aeroelastic system consisting of simple bending and 
torsion modes was used to study aeroelastic behavior at 
transonic speeds. Generalized aerodynamic forces are 
presented for a wide range of Mach number and reduced 
frequency. Aeroelastic characteristics are presented for 
variations in freestream Mach number, mass ratio, and 
bending-torsion frequency ratio. Flutter boundaries are 
presented which have two transonic dips in flutter speed. 
The first dip is the "usual" transonic dip involving a 
bending-dominated flutter mode. The second dip is 
characterized by a single degree-of-freedom torsion 
oscillation. These aeroelastic results are physically 
interpreted and shown to be related to the steady-state shock 
location and changes in generalized aerodynamic forces due to 
freestream Mach number. 
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generalized aerodynamic force resulting from 

pressure induced by mode j acting through 

displacements of mode i 

full-span aspect ratio 

wing root semichord, Cr/2 

wing local chord 

wing reference chord (root chord) 
pressure coefficient 
critical pressure coefficient 
structural mode shape 
reduced frequency, a>Cr/2U 
freestream Mach number 
pulse motion 
pulse amplitude 

nondimensional dynamic pressure, 
(U/(bcD a Vp))2 

nondimensional dynamic pressure at flutter 

Laplace transform variable, a ♦ ia> 

time, nondimensionalized by freestream speed 

and wing reference chord 

nondimensional time at the pulse center 

wing taper ratio 

freestream speed 


w 

constant related to the pulse width 

z(x,y,t) 

wing deformation at time t 

ao 

mean angle of attack 


fractional semispan 

Ac/2 

midchord sweep angle 


ratio of wing mass to mass of air in the 
truncated cone that encloses the wing 

a 

damping coefficient of Laplace transform 
variable 

CO 

angular frequency 

<0h 

natural frequency of bending mode 


natural frequency of torsion mode 


Introduction 

Prediction of transonic unsteady aerodynamics for 
aeroelastic analysis has been aided by the development of 
finite-difference computer codes. An example of one of the 
most fully developed finite-difference codes for aeroelastic 
analysis is the CAP-TSDl code recently developed at the 
NASA Langley Research Center. The name CAP-TSD is an 
acronym for Computational Aoroelasticity Rrogram - 
I/ansonic S.mall Qjsturbance. The code permits the 
calculation of steady and unsteady flows about complete 
aircraft configurations for aeroelastic applications in the 
transonic speed range. Steady and unsteady pressure 
comparisons were presented in Refs. 1 to 3 for numerous 
cases which demonstrated the geometrical applicability of 
CAP-TSD. These calculated results were generally in good 
agreement with available experimental pressure data which 
validated CAP-TSD for multiple component applications with 
mutual aerodynamic interference effects. Preliminary 
flutter applications of CAP-TSD were presented in Ref. 4 for 
a simple well-defined wing case. The calculated flutter 
boundaries compared well with the experimental data which 
gives confidence in CAP-TSD for aeroelastic prediction. 
More recently, the code was modified to include entropy and 
vorticity effects to treat cases with relatively strong shock 
waves. 5 Without these effects, the isentropic formulation 
typically overpredicts the shock strength and locates the 
shock too far aft in comparison with experiment. The 
modified code includes the entropy and vorticity effects while 
retaining the relative simplicity and cost efficiency of the 
small-disturbance potential formulation. As shown in 
Ref. 5, the results obtained by including these effects were 
in very good agreement with parallel Euler calculations for 
attached-flow strong-shock cases. 

In Refs. 1 to 5 the CAP-TSD code was shown to be 
accurate and efficient for transonic unsteady aerodynamic 
and aeroelastic applications. It is therefore timely to utilize 
such a computer code to study the generalized aerodynamic 
forces (GAFs) and aeroelastic stability characteristics of 
wings in transonic flow. The purpose of the paper is 
therefore to report on the results of a transonic unsteady 



aerodynamic and aeroeiastic analysis of a generic transport- 
type wing. The primary objectives of the research were to: 
(1) calculate transonic GAFs and aeroeiastic characteristics 
for a wide range of parameter variations; and (2) correlate 
and physically interpret the resulting pressure 
distributions, GAFs, and aeroeiastic stability characteristics 
including flutter. 

in this study, a binary aeroeiastic system consisting of 
simple bending and torsion modes was used to investigate 
aeroeiastic behavior at transonic speeds. Results were 
obtained for incremental changes in freestream Mach 
number. Effects of mass ratio and bending-torsion 
frequency ratio on aeroeiastic characteristics were also 
considered. The present work is analogous to the numerous 
aeroeiastic studies performed for two-dimensional airfoils 
by many researchers including Isogai fi Yang et al., 7 and 
Edwards, et al.8 In Refs. 6 to 8, aeroeiastic analyses were 
performed for a large range of freestream conditions and 
structural parameter variations for numerous airfoil 
sections. These studies contributed significantly to the 
understanding of aeroeiastic behavior at transonic Mach 
numbers. The present work may similarly give physical 
insight into the flow mechanisms which control transonic 
aeroeiastic behavior of three-dimensional wings. 

The motivations for the study were to gain a better 
understanding of transonic unsteady flows and aeroeiastic 
behavior, and the need to develop efficient computational 
methods to predict and investigate aeroeiastic phenomena. 
For example, the transonic GAFs in the present study were 
obtained by extending the pulse transfer-function analysis of 
Refs. 9 and 10 to treat flexible modes of motion. In this 
anafysis, the GAFs are obtained indirectly from the response 
due to a smoothly varying exponentially shaped pulse. Use of 
the pulse transfer-function analysis gives considerable 
detail in the frequency domain with a significant reduction in 
cost over the alternative method of forced harmonic 
oscillation. The paper presents a description of the pulse 
analysis including results and comparisons which assess the 
capability. Transonic aeroeiastic stability analyses are 
performed with the GAFs from the pulse analysis, using 
state-space aeroeiastic modeling such as that reported in 
Refs. 11 to 14. The stability analysis of the present study is 
an extension of the state-space model of Refs. 13 and 14 for 
binary aeroeiastic anafysis of a three-dimensional wing. A 
brief description of the procedure is given in a subsequent 
section. 

A secondary objective of the present work was to 
determine whether the analysis techniques are capable of 
predicting more challenging aeroeiastic problems of 
significant current interest. For example, the results to be 
presented show a second, more critical (lower) dip in 
flutter speed at transonic Mach numbers slightly higher than 
the first dip. A flutter boundary with two transonic dips has 
indeed been found experimentally for a transport-type 
supercritical-wing flutter model tested at the NLR in the 
Netherlands. The first dip was identified in Ref. 15 as the 
"usual* transonic dip which involves a bending-dominated 
flutter mode, for the wing in predominantly attached flow. 
The second dip, however, occurred in a range of Mach 
number where the flow was separated and was identified as a 
torsional-buzz type of flutter, since the flutter mode was 
almost identical with the torsion vibration mode shape. The 
phenomenon was attributed in Refs. 15 and 16 to an 
aerodynamic resonance which occurred in a narrow range of 
Mach number at higher angles of attack for this wing. It has 
been suggested by Mabey, 17 that if this aerodynamic 
resonance occurs at frequencies close to wing torsional 
frequencies, the single degree-of-freedom torsional buzz 
may result. A similar torsional buzz phenomenon has also 
been reported by Moss and Pierce 1 3 where the primary 
torsion mode of a series of wings tended to be strongly 
excited under separated flow conditions. This excitation led 


to a sustained high-amplitude response or limit-cycle 
torsional buzz, in this case, a strong shock wave occurred on 
the upper surface of the wing, which was approximately 
parallel to the torsion-mode node line. Also, the shock was 
of fairly uniform strength across the span in the outboard 
region of the wing. The movement of the shock with angle of 
attack was such that the slope of the local pitching moment 
about the node line was negative rather than positive. As 
discussed in Refs. 17 and 18, these features seem to be 
necessary conditions for the occurrence of torsional buzz. 

Although the current work is limited to inviscid attached 
flow, the results to be presented exhibit similar flutter 
behavior as observed in the experiments reported in 
Refs. 15 to 18 and described above. Certainly when strong 
shocks are present in the inviscid calculations or when the 
physical situation is known to involve separated flow, the 
present analysis needs to be extended to account for 
unmodeled effects. Nevertheless, the detailed comparisons 
and physical interpretations presented may give an increased 
understanding of unsteady transonic flows and insight into 
the fundamental aeroeiastic mechanisms responsible for 
flutter at transonic speeds. 


Computational Proc«dur»» 

In this section, the computational procedures used in the 
present study are described including the CAP-TSD code, the 
pulse transfer-function analysis, the structural mode 
shapes, and the aeroeiastic stability analysis. 


CAP’TS B.CPdft 

The recently developed CAP-TSD 1 finite-difference code 
was used to calculate the steady and unsteady aerodynamics 
presented in the paper. The CAP-TSD code is an unsteady 
transonic small-disturbance (TSD) code developed for 
transonic aeroeiastic analysis of complete aircraft 
configurations. The code uses a time-accurate approximate 
factorization (AF) algorithm developed by Batina 1 5 for 
efficient solution of the unsteady TSD equation. The AF 
algorithm consists of a Newton linearization procedure 
coupled with an internal iteration technique. The CAP-TSD 
code is capable of treating combinations of lifting surfaces 
and bodies, and includes the following algorithm features: 
(1) Engquist-Osher monotone differencing, (2) 
nonreflecting far field boundary conditions, (3) second- 
order accurate spatial differencing in supersonic regions of 
the flow, and (4) entropy and vorticity effects to treat cases 
with strong shock waves. Further details of the algorithm 
development and solution procedures are reported in Refs. 3, 
5, and 19. A detailed description of CAP-TSD may be found 
in Ref. 1. 


Eulia Transfer-Function Analysis 

Generalized aerodynamic forces which are typically used 
in aeroeiastic analyses may be obtained by calculating 
several cycles of a harmonically forced oscillation with the 
determination of the forces based upon the last cycle of 
oscillation. The method of harmonic oscillation requires one 
flow field calculation for each value of reduced frequency of 
interest. By contrast, GAFs may be determined for a wide 
range of reduced frequency in a single flow field calculation 
by the pulse transfer-function analysis.^* 1 *) In the pulse 
analysis, the GAFs are computed indirectly from the 
response of the flow field due to a smoothly varying 
exponentially shaped pulse. Results computed using the 
pulse analysis for a pitching flat plate airfoil were first 
presented in Ref. 9. These results were in good agreement 
with parallel linear theory calculations which validated the 
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accuracy of the analysis. Applications to transonic airfoil 
cases were also in good agreement with GAFs computed using 
the harmonic method. Calculations for a flat plate wing 
undergoing a simple rigid pitching motion were presented in 
Ref. 10. These three-dimensional results were also shown 
to agree well with linear theory for reduced frequencies less 
than unity. Therefore, because of the computational 
efficiency of the pulse transfer-function analysis, the 
capability was extended in the present study to treat general 
flexible modes such as the natural vibration modes of a wing. 
The pulse is expressed as 

*W(1- t c ) 2 

q(t) = q o e (1) 


where q 0 is the pulse amplitude, w is a constant related to 
the width of the pulse, and tc is the time at the center of the 
pulse. The deformation of any point on the wing at time t, is 
determined by the product of the pulse motion and the mode 
shape. The deformation z(x,y,t) is then given by 

w(t - g 2 

z(x,y,t) = q o c f(x,y) ( 2 ) 


where f(x,y) is the mode shape, and x and y are 
nondimensional Cartesian coordinates in the streamwise and 
spanwise directions, respectively. A small pulse is 
prescribed in a given mode of motion and the aerodynamic 
transients are calculated. The aerodynamic transients are 
then used to obtain the GAFs in the frequency domain by a 
transfer-function analysis. In this analysis, a fast Fourier 
transform (FFT) of the aerodynamic transients is divided by 
the FFT of the corresponding pulse motion to calculate the 
GAFs. The GAFs, Ajj, are defined as the integrated forces 
resulting from pressure induced by mode j acting through 
the displacements of mode i. The transform assumes that the 
system is linear which is a valid approximation for most 
transonic cases, since experience has shown that the 
response for harmonic or aeroelastic motion is, in general, 
locally linear for small amplitudes of osciilation.9’ 14 


Structural Mud* Stiattsa 

in this investigation, simple polynomial equations were 
assumed to describe the vibration mode shapes. The 
following equations simulate first bending and first torsion 
mode shapes of a wing, respectively: 

f B = y [(x - j) sin A c/2 + y cos A c/2 ) ( 3 ) 


£r = y [(* - y) cos \ n - y sin A c/2 ] ( 4 ) 


where \ c n is the midchord sweep angle of the wing. The 
equation describing the mode shape due to bending (Eq. (3)) 
assumes that a node line is perpendicular to the wing 
midchord line at the root. The equation describing the mode 
shape due to torsion (Eq. (4)) assumes that a node line 
coincides with the wing midchord line. The amplitude of the 
pulse motion is selected to give a maximum tip deflection due 
to bending of approximately four percent of the semispan and 
a maximum angle of attack at the tip due to torsion of one 
degree. 


Afttflftlaallfi. Stability Analyse 

Constant sectional properties from root to tip are 
assumed for the aeroelastic analyses. This simplified 
modeling is consistent with the assumed mode shapes of Eqs. 
(3) and (4), which in turn assumes that the aeroelastic 
motion of the wing can be described by a linear combination 
of the fundamental bending and torsion modes. The 
calculations were performed using state-space aeroelastic 
modeling such as that reported in Refs. 11 to 14. Edwards, 
et al. 1 1 used a state-space model employing Fade' 
approximants to model the unsteady airloads and 
demonstrated good agreement with a time-marching 
technique for a linearized case. Bland and Edwards 1 2 
demonstrated that such locally linear procedures may be 
used with airloads derived from a transonic small- 
disturbance code for airfoil applications. Batina and 
Yang 1 3. 1 4 used a similar procedure to study transonic 
aeroelastic stability and response behavior of airfoils with 
active controls. The stability analysis of the present study is 
an extension of the state-space model of Refs. 13 and 14 for 
binary aeroelastic analysis of a three-dimensional wing. 
The resulting locally linear stability model provides a 
relatively inexpensive determination of aeroelastic 
stability, while retaining the nonlinear properties of the 
transonic mean flow. The model is derived by assuming a 
linear superposition of airloads due to bending and torsion 
motions. The required airloads are approximated by curve- 
fitting the CAP-TSD GAFs with a Pade' approximating 
function.^ The function may then be rewritten as a set of 
ordinary differential equations which, when coupled to the 
equations of motion and Laplace transformed, leads to a 
linear first-order matrix equation 

— {2)=[A]{2) (5) 

0) 

a 

where {z} contains the displacements, velocities, and 
augmented states and (A] is a real matrix of constant 
elements. Equation (5) is solved using linear eigenvalue 
solution techniques for specified values of the 
nondimensional dynamic pressure Q, defined as 
(U/(bo> a >/jI))2. The resulting eigenvalues are plotted in a 
dynamic pressure "root- locus” type format in the complex 
s-plane. 


BaiuHa and Discussion 

Results are first presented from the pulse transfer- 
function analysis to assess the accuracy and efficiency of the 
capability by making comparisons with other calculations. 
The pulse analysis was then used in a parameter study to 
investigate the effects of freestream Mach number on GAFs at 
transonic speeds. Results are also presented from the 
aeroelastic stability analysis, obtained using the GAFs from 
the pulse analysis. 


Cflflfiqu ration.- Analyzed 

The wing selected for the present study is the Royal 
Aircraft Establishment (RAE) research wing *A" that was 
tested in the RAE 8ft. x 6ft. transonic wind tunnel. 2 0 The 
wing is a structurally rigid model that has a simple 
planform without dihedral or twist and a symmetric RAE 
101 airfoil section with a maximum thickness- to-chord 
ratio of 9%. As shown in Fig. 1, the wing has a midchord 
sweep angle of 30°, a taper ratio of one-third, and a full- 
span aspect ratio of six. The polynomial mode shapes that 
were assumed for the wing are shown in Fig. 2. The RAE 
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wing A is a standard configuration selected by the AGARD 
Fluid Dynamics Panel for steady flow calculations. Although 
the wing is a rigid model and no flutter data is available, the 
wing was selected for the present study because of the 
simplicity of the planform and the availability of 
experimental steady pressure data. 20 The experimental 
steady pressure data that is reported in Ref. 20, however, 
was measured with the wing mounted on an axisymmetric 
fuselage. Therefore, a simple fuselage was modeled with 
CAP-TSD, along with the wing, to allow for direct 
comparison with the steady pressure data. All of the other 
calculations to be presented were performed without the 
fuselage. 

Calculations were performed for thirteen values of 
freestream Mach number including M - 0.8, 0.825, 0.85, 
0.875, 0.9, 0.91. 0.92, 0.93, 0.94, 0.95, 0.96, 0.97, and 
0.98. Results are presented mainly for a subset of these 
cases, namely at M * 0.8, 0.9, 0.92, 0.94, 0.96, and 0.98. 
For M 2 0.94, relatively strong shock waves are present in 
the flow. For these cases, the isentropic calculation 
typically overpredicts the shock strength and locates the 



RAE wing "A" 


a c/2 -30°. AR -6.0, TR - 1/3 
Fig. 1 RAE wing "A" geometry. 


Mode 1 : Bending 



Mode 2: Torsion 



Fig. 2 Structural mode shapes. 


shock too far aft in comparison with the physical situation. 
Therefore, the calculations were repeated for M 2 0.94 by 
including the entropy and vorticity effects to accurately 
treat these cases. Furthermore, calculations were also 
performed at M - 0.9 for a series of wings created by 
varying the midchord sweep angle from A c /2 - 14° to 40°. 
These results were obtained to determine effects of wing 
sweep on transonic GAFs and flutter. These results may be 
found in Ref. 21 . 


Steady EttiiiUft Baiulli 

Steady pressure results are presented first to assess the 
basic character of the transonic flow. Since the unsteady 
flow field and hence generalized aerodynamic forces depend 
strongly on the transonic steady-state flow, the steady 
pressure distributions frequently can give physical insight 
into critical flow mechanisms which control aeroelastic 
phenomena. Of interest are the flow characteristics in the 
wing tip region such as the shock strength and shock location 
relative to the node line of the torsion mode. 


CAP-TSD 



(a) i\ - 0.40. 


CAP-TSD 



(b) q - 0.85. 

Fig. 3 Comparison between CAP-TSD and experimental 
steady pressure distributions on the RAE wing 
"A" at M - 0.90 and ao - 0°. 
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M = 0.80 0.90 0.92 0.94 0.96 0.98 



Experimental Data Comparison . - To assess the 
accuracy of the CAP-TSD solution for this configuration, 
results were obtained for the RAE wing A and fuselage at 
M • 0.9 and ao - 0° for comparison with the experimental 
data. Steady pressure distributions at two span stations 
along the wing are presented in Fig. 3. Results for f\ - 0.40 
are shown in Fig. 3(a); results for f\ - 0.85 are shown in 
Fig. 3(b). The calculated pressures are identical along the 
upper and lower wing surfaces, since the airfoil section is 
symmetric and the wing is at 0° angle of attack. These 
pressure distributions indicate the presence of embedded 
supersonic regions along the upper and lower surfaces. A 
mild shock wave also occurs in the inboard region of the wing 
(not shown). In general, the CAP-TSD pressures agree well 
with the experimental data. 


Eff.ftc.ta of Freealreim. Mach . Number - steady 

pressure distributions at six different Mach numbers and 
ao - 0° are presented in Fig. 4. The pressures for M > 0.94 
were computed including the entropy and vorticity 
corrections. At M * 0.8, the calculated pressures indicate 
that the flow is entirely subcritical. At M - 0.9, there is a 


TSD ♦ entropy ♦ vorticity 



Fig. 5 Steady pressure distributions at n - 0.9 for a 
range of Mach number at ao - 0°. 


mild shock wave in the inboard region of the wing that 
weakens in the outboard direction toward the wing tip. When 
the freestream Mach number was increased to M « 0.92 and 
0.94, a moderately strong shock wave forms across the span 
along the upper and lower wing surfaces. In the outboard 
region of the wing, the shock is approximately parallel to the 
torsion-mode node line at the local midchord. Here the shock 
is also of fairly uniform strength in the spanwise direction. 
These features of the shock location and strength are similar 
to those reported in Refs. 17 and 18, which were considered 
to be necessary (but not sufficient) conditions for torsional 
buzz to occur. When the Mach number was further increased 
to M - 0.96 and 0.98, the shock moves aft towards the 
trailing edge. Since the steady pressure distributions show 
relativeiy small changes in the range 0.8 s M s 0.9, it is 
anticipated that only small changes will occur in the GAFs. 
Between M - 0.92 and 0.98, much larger changes in the 
GAFs are expected because of the larger differences in the 
steady flows, especially in the wing tip region. 

Figure 5 shows more detailed comparisons of the steady 
pressure distributions at fj - 0.9 for the six Mach 
numbers ranging from M - 0.8 to 0.98. These comparisons 
of pressures clearly show the transonic steady flow 
characteristics in the wing tip region. For M 2 0.94, 
results corresponding to the calculations with and without 
entropy and vorticity effects are shown. Without entropy 
and vorticity effects, the shock is stronger and moves aft 
much faster with increasing Mach number. At M - 0.96 and 
0.98, for example, the isentropic shock is located 
approximately 20% chord downstream of the nonisentropic 
shock. Also, at M » 0.98 the isentropic calculation predicts 
that the shock wave is at the trailing edge. Furthermore, the 
pressure distributions of Fig. 5 show that the shock wave 
crosses the torsion-mode node line with increasing Mach 
number near M - 0.94. This fact is expected to play a 
significant role in the aeroelastic behavior of the wing. 


Genwalltad Aerodynamic Force Ra«ult« 

Generalized aerodynamic forces are presented beginning 
with an assessment of the accuracy and efficiency of the 
pulse transfer-function analysis relative to harmonic 
oscillation calculations. Further results are then presented 
from the calculations which include the nonisentropic 
corrections to demonstrate the effects of freestream Mach 
number on GAFs. All of the results are presented in the form 
of real and imaginary components of the GAFs, Ay, as 
functions of the reduced frequency k. In the present study, 
bending and torsion are defined as modes 1 and 2, 
respectively. 


5 




LlOfttr Theory CttmgarllQn- - To assess the accuracy 
of the pulse transfer-function analysis to treat flexible mode 
shapes, results were obtained for the wing modeled as a flat 
plate at M - 0.9 and ao - 0°. These results are compared 
with linear theory results computed using RHOIV, an 
assumed pressure mode kernel function method for simple 
harmonic motion 22,23 The RHOIV generalized aerodynamic 
forces were obtained for eleven values of reduced frequency 
ranging from k « 0.0 to 1.0 in increments of 0.1. 
Comparisons of results from the pulse analysis and RHOIV 
are given in Fig. 6. The results shown are for the bending 
mode GAF due to torsion motion, Ai2> which is similar to 
wing unsteady lift-curve slope. The two sets of results are 
in good agreement, for the wide range of k shown, which 
validates the pulse analysis for application to flexible mode 
shapes. Results for the other three QAFs show similar good 
agreement between CAP-TSD and RHOIV and therefore are not 

shown. 21 


Transonic Comparison . - To assess the accuracy of 
the pulse analysis for transonic applications, GAFs were 
obtained for the wing at M - 0.9 and oo - 0°. The accuracy 
of these forces is determined by making comparisons with 
similar results computed using the harmonic method. 
Results from the harmonic method were obtained for eleven 
values of reduced frequency ranging from k - 0.0 to 1.0 in 
increments of 0.1. Comparisons of the A12 GAF calculated 
by the pulse analysis and the harmonic method are shown in 
Fig. 7, The two sets of results are in very good agreement 
which validates the pulse transfer-function analysis for 
application to transonic cases. Results for the other three 
GAFs show similar good agreement between the pulse 
analysis and the harmonic method and therefore are not 

shown. 21 

With respect to computational efficiency, the harmonic 
method required one flow field calculation for each value of 
reduced frequency for each mode of motion. This is in 
contrast to the pulse analysis which determined the GAFs for 
the entire range of frequency in a single flow field 



Fig. 6 Comparison between generalized aerodynamic 
forces calculated by CAP-TSD (linear) and 
RHOIV for the wing modeled as a flat plate at 
M « 0.9 and ao - 0°. 


calculation for each mode. Since eleven values of reduced 
frequency were selected for calculation by the harmonic 
method, the computational cost of obtaining the pulse results 
was consequently about an order of magnitude less than the 
cost of the harmonic method results. For the case considered 
here, the pulse analysis is therefore a relatively 
inexpensive method of determining transonic GAFs which in 
turn makes it economically feasible to conduct the parameter 
variation study reported below. 


Effects of Fleet tr earn _ M a c h Number . - To 

determine the effects of freestream Mach number on the 
GAFs, results are presented for the wing at M - 0.8, 0.9, 
0.92, 0.94, 0.96. and 0.98 at zero mean angle of attack, 
corresponding to the steady pressure distributions presented 
in Figs. 4 and 5. Comparisons of GAFs at these Mach 
numbers are presented in Figs. 8(a) to 8(d). 

For the GAFs resulting from the pressure induced by 
bending (An and A21 shown in Figs. 8(a) and 8(b), 
respectively), the real part of A21 shows the most 
significant change due to Mach number. With increasing 
Mach number, the real part of A21 decreases in magnitude 
for 0.92 s M < 0.96. This decrease is attributed to an aft 
shift in the aerodynamic center which begins to occur after 
M - 0.92 due to the formation of the shock wave in the tip 
region. As shown in the steady pressure distributions of 
Figs. 4 and 5, there is a mild shock wave at M - 0.92 which 
grows in strength when the Mach number is increased to 
M - 0.94. Since A21 is similar to a wing moment coefficient 
and the torsion-mode node line is located at the local 
midchord, an aft shift in the aerodynamic center produces a 
decrease in magnitude of the real part of A21. Significant 
changes in aeroeiastic behavior may occur for M > 0.92 
because of this large decrease in magnitude of A21 . 

For the GAFs resulting from the pressure induced by 
torsion (Ai 2 and A22 shown in Figs. 8(c) and 8(d), 
respectively), increasing the Mach number produced large 
changes in the forces especially in the low k range. These 



Reduced frequency k 


Fig. 7 Comparison between CAP-TSD generalized 
aerodynamic forces calculated by pulse analysis 
and harmonic analysts for the wing at M • 0.9 
and ao • 0*. 



M 


M 




0.80 

0.90 

0.92 

0.94 

0.96 

0.98 


(a) GAF resulting from pressures due to 
bending weighted by bending displacements. 


(c) GAF resulting from pressures due to 
torsion weighted by bending displacements. 


0.80 

0.90 

0.92 

0.94 

0.96 

0.98 



0.80 

0.90 

0.92 

0.94 

0.96 



* 


V 


(b) GAF resulting from pressures due to 
bending weighted by torsion displacements. 


(d) GAF resulting from pressures due to 
torsion weighted by torsion displacements. 


Fig. 8 Comparison of generalized aerodynamic forces 
(GAFs) for a range of Mach number at ao • 0°. 
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changes in GAFs are of primary interest because they occur 
in the range of reduced frequency where flutter typically 
occurs (0.1 < k < 0.4). With increased Mach number from 
M - 0.8 to 0.92, for example, both the real and the 
imaginary parts of At 2 increase in magnitude. The At 2 
force is similar to wing lift-curve slope which increases as 
Mach number increases. This increases both the magnitude 
and the phase lag of A 12 which in turn controls the 
aeroelastic behavior for 0.8 £ M £ 0.92. As discussed by 
IsogaiG and later by Zwann 24 it is the increased phase lag of 
A 1 2 which is the flow mechanism that produces the 
transonic dip in flutter speed. Also with increased Mach 
number, the real part of A 22 generally decreases which is 
again attributable to an aft shift in aerodynamic center due to 
shock waves. Furthermore, for the higher Mach numbers in 
the range 0.94 s M £ 0.98, the imaginary part of A 22 
changes sign for low reduced frequencies. This change in 
sign indicates the possibility of single degree-of-freedom 
flutter in torsion. 


Aeroelastic Stability and Flutter Result* 

Aeroelastic stability and flutter results are presented 
next to investigate aeroelastic characteristics at transonic 
speeds. Calculations were performed for several values of 
freestream Mach number, mass ratio, and for a>h/a>a - 0.1, 
0.3, and 0.5. Root-loci plots for p * 100 and o)h/co a - 0.3 
are presented in Fig. 9 for the six Mach numbers 
corresponding to the steady pressure distributions of Figs. 4 
and 5 and the GAFs of Figs. 8(a) to 8(d). Stability results 
for M - 0.8, 0.9, 0.92, 0.94, 0.96, and 0.98 are shown in 
Figs. 9(a) through 9(f). respectively, computed including 
the nonisentropic corrections. The stability plots show the 
migration of bending and torsion root-lod as Mach number 
is increased. The increment in dynamic pressure shown in 
the figures is aQ - 0.2. At M - 0.8, the curves indicate a 
classical bending-torsion type of aeroelastic behavior. Here 
the torsion branch moves to the left in the stable left-half 





(a) M - 0.80. 


(b) M - 0.90. 


(C) M - 0.92. 



(d) M - 0.94. 


(e) M - 0.96. 


(0 M - 0.98. 


Fig. 9 Nondimensional dynamic pressure root-loci for a range of 
Mach number at ao - 0°, p • 100, coh/coa * 0.3, and AQ * 0.2. 
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plane, with increasing dynamic pressure, while the 
bending-dominated branch becomes the flutter mode. With 
increasing Mach number in the range 0.8 £ M s 0.92, the 
flutter dynamic pressure decreases and flutter occurs at a 
lower frequency. Between M - 0.92 and 0.94, the rooMoci 
coalesce in frequency and then switch, in that the flutter 
crossings at M - 0.94, 0.96, and 0.98 originate from the 
wind-off torsion mode rather than the bending mode. A rapid 
) increase in the flutter frequency thus occurs in the range 

0.92 s M s 0.96, from a value close to the bending natural 
frequency at M * 0.92 to a value close to the torsion natural 
frequency at M « 0.96. These aeroelastic characteristics are 
further explained by considering the corresponding flutter 
boundaries which are given in Fig. 10, and correlating these 
results with the steady pressures of Figs. 4 and 5, and the 
GAFs of Figs. 8(a) to 8(d). Figure 10(a) shows the flutter 



M 


speed index (defined as U/(bco a ^)) versus Mach number 
and Fig. 10(b) shows the nondimensional flutter frequency 
(defined as o)/a> a ) versus Mach number. In these figures, 
the points along the p - 100 curves which are labeled (a) 
through (f), correspond to the root-loci of Figs. 9(a) to 
9(f). respectively. Both isentropic as well as nonisentropic 
curves are plotted. Flutter results are also shown in Fig. 10 
for n - 20 to indicate effects due to mass ratio. 

The flutter boundaries of Fig. 10(a) show two dips in 
flutter speed index, the second of which is deeper than the 
first. These transonic dips in flutter speed are deeper at the 
higher mass ratio of n - 100. More specifically, as Mach 
number increases in the range 0.8 s M £ 0.92, the flutter 
speed decreases as expected. This decrease is due to the 



(a) flutter speed index versus Mach number. 



M 


(b) nondimensional flutter frequency versus 
Mach number. 

Fig. 10 Effects of mass ratio on flutter characteristics at 
ao * 0°. 


(a) flutter speed index versus Mach number. 



M 


(b) nondimensional flutter frequency versus 
Mach number. 

Fig. 11 Effects of bending-torsion frequency ratio on 
flutter characteristics at ao - 0°. 
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increased phase lag of the A 12 generalized aerodynamic force 
which occurs with increasing Mach number (Fig. 8(c)). 
The bottom of the first dip occurs at M - 0.91 - 0.92 for 
both pi - 20 and 100. Here the shocks are relatively weak, 
and the frequency curves of Fig. 10(b) suggest that the 
flutter mode is predominately bending since the flutter 
frequency is close to the bending mode frequency. Inspection 
of the corresponding flutter eigenvectors confirmed this. At 
M - 0.92 and p - 100, for example, the ratio of the 
magnitudes of bending and torsion in the flutter mode 
eigenvector is 8:1. With a further small increase in Mach 
number, there is a rapid increase in flutter speed 
(Fig. 10(a)) which defines the "other side" of the first dip. 
Physically this occurs because there is a significant decrease 
in the magnitude of the A 21 and A 22 GAFs for 
0.92 < M < 0.94 (Figs. 8(b) and 8(d)) due to the aft shift in 
aerodynamic center. Recalling that these GAFs are tike 
moment coefficients, they decrease in magnitude when the 
aerodynamic center shifts aft toward the torsion-mode node 
line. Associated with this change in GAFs is a significant 
decrease in the magnitude of the imaginary part of A 22 
which indicates a loss of aerodynamic damping in torsion. 
The decrease in torsional damping causes the flutter mode to 
rapidly transition from a bending-dominated instability at 
M - 0.92 to a torsion-dominated instability at M - 0.94. 
This transition is shown by the rapid increase in flutter 
frequency from M - 0.92 to 0.94 shown in Fig. 10(b). It 
was also shown in Figs. 9(c) and 9(d), by the switch in 
root-loci, which results in a much higher frequency at the 
flutter crossing. Near M « 0.94, the shock, which is 
approximately parallel to the torsion-mode node line in the 
outboard region of the wing, crosses the node line with 
increasing Mach number. This has a destabilizing effect with 
respect to flutter since the damping in torsion is actually 
negative for low values of k as shown in Fig. 8(d). The 
flutter speeds in the first part of the second dip are lower in 
the isentropic calculations since the shock wave is stronger 
and located further aft of the torsion-mode node line, as 
previously shown in Fig. 5. The further loss of torsional 
damping at the higher Mach numbers produces a flutter mode 
that is almost identical to the torsion mode shape, as 
suggested by the flutter frequency values for M > 0.95 shown 
in Fig. 10(b). The ratio of the magnitudes of torsion and 
bending in the flutter eigenvector at M - 0.97 and p - 100, 
for example, is 80:1 . Therefore the instability in the second 
dip is a single degree-of-freedom flutter similar to the 
torsional-buzz phenomenon reported in Refs. 15 to 18. 
With a further increase in Mach number, the shock wave in 
the isentropic calculations moves aft to the trailing edge near 
M - 0.97, which has a stabilizing influence on the 
aeroelastic system since the flow about the wing is mostly 
supersonic. With the shock at the trailing edge, the flutter 
speed increases rapidly thus defining the "other side" of the 
second flutter dip. in the more-accurate nonisentropic 
calculations, however, the shock wave does not reach the 
trailing edge for the range of freestream Mach number 
considered. Hence, the flutter characteristics are 
considerably different and the importance of including 
entropy and vorticity effects is evident. 

Effects of bending-torsion frequency ratio on flutter 
characteristics are shown in Fig. 11. Flutter speed index 
and nondimensional flutter frequency as functions of 
freestream Mach number are shown in Figs. 11(a) and 
11(b), respectively. Calculations were performed for 
coh/wa - 0.1, 0.3, and 0.5 with a mass ratio of p - 100. In 
Fig. 11, the curves for ©h/© a - 0.3 and p - 100 are the 
same as those shown in Fig. 10. The results indicate that the 
first dip in flutter speed is significantly influenced by 
o)h/(o a - With a)h/a>a - 0-5, for example, the dip is very 
shallow in comparison to the dip with ©h/©a ■ 0.1. By 
decreasing ©h/©a from 0.5 to 0.1, the flutter speed 
decreases by 56%. The bottom of the dip also shifts to a 


slightly higher Mach number. The flutter frequency curves 
of Fig. 11(b) show the same trends for the three frequency 
ratios, each curve being shifted toward the respective 
bending mode frequency for M < 0.93. The flutter mode at 
the bottom of the first flutter dip, however, is somewhat 
different for the three cases. Inspection of the flutter 
eigenvectors revealed that the ratios of magnitudes of 
bending and torsion for ©h/©a - 0.1, 0.3, and 0.5 are 70:1, 
8:1, and 5:1, respectively. This indicates that the flutter 
mode at the bottom of the dip for ©h/©a - 0.1 is a single 
degree-of-freedom bending instability. The results further 
indicate that the second dip in flutter speed is independent of 
the bending-torsion frequency ratio which is not unexpected 
since here the flutter mode is a single degree-of-freedom 
torsional instability. 


Concludlno_RemarK« 

Transonic aeroelastic stability analyses and flutter 
calculations were presented for a generic transport-type 
wing based on the use of the CAP-TSD (Computational 
Eeroelasticity Erogram - Iransonic Email disturbance) 
finite-difference code. In this study, a binary aeroelastic 
system consisting of simple bending and torsion modes was 
used to study aeroelastic behavior at transonic speeds. 
Transonic generalized aerodynamic forces (GAFs) and 
aeroelastic characteristics were computed for a wide range 
of freestream Mach number. The CAP-TSD calculations were 
performed by including entropy and vorticity effects to more 
accurately treat cases with strong shock waves. 

An efficient method of calculating transonic generalized 
aerodynamic forces was extended to treat general flexible 
modes of motion of a wing. The method presented was the 
pulse transfer-function analysis which determines the 
unsteady forces indirectly from the response due to a 
smoothly varying exponentially shaped pulse. Comparisons 
of calculated GAFs for both linear and nonlinear cases showed 
good agreement with alternative calculations which verified 
the pulse analysis for application to flexible modes. 
Furthermore, the computational resources required by the 
pulse analysis were an order of magnitude less than those 
required by the harmonic method. 

The pulse analysis was used in a parameter study to 
investigate the effects of freestream Mach number on GAFs at 
transonic speeds. The GAFs from the pulse analysis were 
utilized in aeroelastic stability and flutter calculations. 
Flutter boundaries were computed which had two transonic 
dips in flutter speed. The first dip was the "usual" transonic 
dip which involved a bending-dominated flutter mode. The 
second dip was characterized by a single degree-of-freedom 
torsion oscillation. This single degree-of-freedom flutter 
occurred at the higher Mach numbers considered due to a loss 
of aerodynamic damping in torsion when the steady shock 
location was downstream of the torsion-mode node line. 
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